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vertical  motion  of  pollutants  in  the  stratosphere.  One  method  to  estimate 
vertical  transport  due  to  the  effects  of  turbulence  is  by  means  of  the  "effective 
diffusivity*"ycoefficient.  To  calculate  this  parameter,  it  is  often  necessary  to 
know  the  value  of  e,  the  rate  of  turbulent  dissipation.  This  parameter,  e,  is 
also  important  to  know  in  the  context  of  the  global  numerical  stratospheric 
simulation  models  now  being  created  for  environmental  assessment  purposes. ? 
One  of  the  most  often  used  methods  to  measure  e employs  turbulence  spectra. 
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20.  Abstract  (Continued) 

These  spectra  as  seen  in  the  HICAT  data,  however,  have  anomalous  properties 
which  must  be  understood  before  we  can  intelligently  use  them  outside  the 
isotropic  range  for  the  purpose  of  measuring  e.  The  HICAT  power  spectra  of 
stratospheric  turbulence  had  -5/3  slopes  on  log-log  graphs  that  extended  up  to 
wavelengths  as  long  as  40,  000  ft.  Standard  turbulence  theory  of  the  inertial 
subrange  would  account  for  this  observation  only  to  wavelengths  of  about  150  ft. 
The  purpose  of  this  report  is  to  extend  existing  theory  to  explain  this  phenom- 
enon. The  results  should  apply  to  all  turbulence  generated  by  the  Kelvin- 
Helmholtz  instability  in  a stratified  fluid  far  from  a boundary.  They  also  may 
explain  certain  well  known  spectral  nonlinearities,  but  the  most  important 
practical  outcome  is  that  a spectral  estimate  of  turbulent  dissipation  rate  must 
be  based  on  itotropic  range  wave  numbers  to  avoid  error  despite  the  -5/3 
slope. 
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Theoretical  Explanation  of  Spectral  Slopes  in 
Stratospheric  Turbulence  Data  and 
Implications  for  Vertical  Transport 


I.  INTRODUCTION 

Stratospheric  turbulence,  like  tropospheric  CAT,  takes  the  form  of  layers 

which  are  broad  in  the  horizontal  direction  (for  example,  25  km)  but  thin  in  the 

vertical  direction  (typically  between  100  m to  1000  m depending  on  how  measured 
12  3 

(Rosenberg  and  Dewan,  Anderson,  and  Barat.  The  power  spectra  of  the  turbu- 
lent velocity  fluctuations  within  such  layers  have  slopes  (on  log-log  plots)  which 
are  roughly  in  the  vicinity  of  -5/3  to  wavelengths  as  long  as  10,000  ft,  or  possibly 
as  long  as  40,  000  ft.  (See  Figure  1.  ) More  specifically,  these  slopes  fall  in  the 
range  of  -1.5  to  -1.7  for  horizontal  and  -1.25  to  -1.4  for  vertical  fluctuations 
(HICAT  data*).  While  one  would  expect  slopes  of  -5/3  for  spectra  in  the  inertial 
sub-range,  one  would  not  ordinarily  expect  to  find  them  at  wavelengths  far  exceed- 
ing the  scale  for  isotropy.  Isotropy  is  found  up  to  100  to  150  ft  in  the  case  of 
stratospheric  turbulence  (HICAT  data*).  The  purpose  of  this  paper  is  to  explain 
the  presence  of  this  anomalous  -5/3  range  of  spectra  in  data  from  stratospheric 

(Received  for  publication  30  September  197R) 

♦References  to  HICAT  data  are  found  in  Crooks,  et  al,  ref.  45.  ) 

1.  Rosenberg,  N.W.,  and  Dewan,  E.M.  (1975)  Stratospheric  Turbulence  and 

Vertical  Effective  Diffusion  Coefficients.  AFCRL-TR-75-0519. 

2.  Anderson,  A.D.  (1957)  Free-Air  turbulence,  J.Meteorol.  14 :477-494. 

3.  Barat,  J.  (1975)  Etude  experimentale  de  la  structure  du  champ  de  turbulence 

dans  la  moyenne  stratosphere,  C.  R.  Acad.  Sc.  Paris  280(Serial  B):691-693. 
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INVERSE  WAVELENGTH  - cycles/ft 


Figure  1.  Typical  HICAT  Turbulence  Spectrum 
Showing  -5/3  Slope  Beyond  the  Isotropic  Range 

1 

turbulence.  As  mentioned  in  the  abstract,  this  is  of  use  in  the  general  vertical 
transport  problem. 

It  is  generally  agreed  that  CAT,  and  presumably  most  stratospheric  turbu- 
lence, is  generated  via  the  Kelvin-Helmholtz  (K-Il)  instability  mechanism 
1 4 5 

(Rosenberg  and  Dewan,  and  Dutton  ' ).  When  the  shear  across  a stratified  layer 
■ is  sufficient  to  cause  the  Richardson  number  to  go  below  0.25,  an  unstable  wave 

j presumably  builds  up  and  eventually  breaks  down  in  such  a way  as  to  generate  a 

layer  of  turbulence.  Excellent  experimental  information  and  photographs  will  be 


4.  Dutton.  J.A.  (1971)  CAT  aviation  and  atmospheric  science,  Rev.Geophys. 

9:613-657. 

5.  Dutton,  J.A.  (1973)  Recent  perspectives  on  turbulence  in  the  free  atmosphere, 

in  N.  K.  Vinnichenko,  N.  Z.  Pinus,  S.  M.  Shmeter,  and  G.  N.  Shur, 
Turbulence  in  the  Free  Atmosphere,  Consultants  Bureau,  New  York. 
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found  in  the  work  of  Thorpe,  Woods,  ' Battan,  and  Ludlam.  The  classic 
theoretical  paper  is  that  of  Rosenhead.  1 1 

The  K-H  mechanism  is  responsible  for  turbulence  to  be  found  below  the  sur- 
face of  the  ocean,  and  it  has  been  said  that  the  resulting  thin  layers  of  turbulence 
which  are  thus  generated  intermittently  in  space  and  time,  are  the  primary  cause 
of  what  vertical  transport  there  is  in  the  interior  of  the  world  ocean  (Woods  and 
Wiley12).  It  is  interesting  to  note  that,  in  this  connection,  some  stratospheric 
and  other  physicists  have  become  aware  that  there  is  a useful  and  close  analogy 

1 O 

between  the  dynamics  of  the  ocean  and  the  stratosphere;  for  example,  see  Thorpe, 
Rosenberg  and  De  van,  1 Woods  and  Wiley,  12  and  Stewart  and  Bolgiano.  14 

An  anomalous  -5/3  slope  (that  is,  one  found  outside  the  inertial  range)  was  not 
seen  for  the  first  time  in  the  HICAT  spectral  data  by  any  means,  but  rather  it  has 
been  repeatedly  seen  in  stratified  media.  For  example,  Pond  et  al1  J reported  the 
results  of  an  experiment  where  the  theoretical  lower  limit  of  the  inertial  range 
was  k 0.  1 cm*1  (k  is  the  wavenumber)  whereas,  in  contrast,  the  k 1 1 spectrum 
provided  a reasonable  fit  to  k 0.005  cm"1!  In  another  case,  Stewart1  ’ made  the 
following  statement  which  is  perhaps  the  most  precise  description  of  the  problem  at 
hand:  "It  is  well  known  that  -5/3  spectra  appear  very  frequently  when  one  observes 
along  a horizontal  path  in  a stratified  fluid.  These  spectra  are  observed  even  when 
it  is  quite  apparent  that  the  turbulence  observed  is  very  far  from  isotropic.  • • • 

This  point  needs  more  carefull  consideration  • • • • than  it  has  yet  received.  " 

0.  Thorpe,  S.  A.  (1973b)  Experiments  on  stability  and  turbulence  in  a stratified 

shear  flow,  J.  Fluid  Alech.  fU(Part  4):731-751. 

7.  Woods,  J.D.  (1969)  On  Richardson's  number  as  a criterion  for  laminar- 

turbulent -laminar  transition  in  the  ocean  and  atmosphere.  Radio  Sci. 

L2:1289-1298. 

8.  Woods,  J.D.  (1968)  Wave-induced  shear  instability  in  the  summer  thermodine, 

J.  Fluid.  Mech.  32:791-800. 

9.  Battan,  L.J.  ( 1 973)  Radar  Observation  of  the  Atmosphere,  l niversitv  of 

Chicago  Press. 


10. 

11. 

12. 

13. 

14. 

15. 

16. 


Ludlam,  F.  If.  (1967)  Characteristics  of  billow  clouds  and  their  relation  t 
clear  air  turbulence.  Quart.  J . Roy.  Met.  Soc.  93:419-43  i. 

Rosenhead,  L.  (1931)  The  formation  of  vortices  from  a surface  of  discontin- 
uity, Proc.  Roy.  Soc.  134:170-192. 

Woods,  J.D. , and  Wiley,  R.  L.  (1972)  Billow  turbulence  and  ocean  micro- 
structure,  Deep  Sea  Research  and  Oceanic  Abst.  19:87-121. 

Thorpe,  S.A.  (1973a)  CAT  in  the  Laboratory,  Weather. 

Stewart,  R.W.  and  Bolgiano,  R.  (1969)  Comments  in:  Clear  air  turbulence 
and  its  detection,  edited  by  Y.H.  Pao  and  A.  Goldberg,  Plenum  Press, 
p.  515-520. 

Pond,  S. , Steward,  R.W. , and  Burling,  R.W.  (1963)  Turbulence  spectra  in 
the  wind  over  waves,  J . Atm.  Sci.  20:319-321. 

Stewart,  R.W.  (1969)  Turbulence  and  waves  in  a stratified  atmosphere. 
Radio  Sci.  4:1269-1278. 
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In  another  case,  Pao  described  this  phenomenon  in  a stratified  turbulent 
flow  as  seen  in  a towing  tank  as  follows:  "However,  the  turbulence  is  too  weak  to 
have  a distinct  f‘D/ J inertial  range.  This  'anomalous  -5/3  subrange,  ' we  believe, 
is  due  to  the  presence  of  the  internal  wave  spectral  peak,  lifting  up  the  spectral 
curve  as  a result.  . . . However,  our  results  (Figure  3)  indicate  that  the  presence 
of  the  f-5/^  frequency  range  does  not  warrant  the  existence  of  an  inertial  subrange 
in  a stably  stratified  fluid."  Figure  3 was  a spectrum  with  -5  3 slope  on  the  log- 

log  plot  showing  the  anomalous  -5/3  subrange. 

18 

Finally,  in  this  connection,  Dutton  wrote  the  following:  "One  of  the  impor- 
tant properties  of  atmospheric  turbulence  is  that  the  —5/3  inertial  subrange  does 
appear  to  exist,  but  for  reasons  unknown  the  -5/3  relation  holds  at  much  smaller 
wave  numbers  than  would  be  expected  from  the  theory." 

From  the  previous  information,  it  is  clear  that  the  anomalous  -5  3 effect  in 
stratified  fluids  is  well  known  and  that  a theoretical  explanation  for  it  would  be  of 
value  at  this  time.  (One  further  mention  of  this  subject  will  be  found  in 
Ellison.  1 

In  order  to  gain  some  appreciation  of  how  far  the  HICAT  spectra  differ  from 
the  usual  predictions  of  the  theory  for  the  inertial  subrange,  we  now  make  a com- 
parison. In  the  case  of  shear  generated  turbulence  in  the  absence  of  buoyancy,  the 
maximum  wavelength  for  isotropy  is  derived  by  comparing  the  mean  rate  of  strain 
to  the  rate  of  strain  of  the  turbulence.  The  latter  is,  from  a similitude  argument, 
given  by  (k3  b(k))1  ^ 2 /2t  where  p(k)  is  the  three-dimensional  scalar  wavenumber 
spectrum  of  the  velocity  fluctuations.  Lettings  be  the  scale  of  the  mean  velocity, 
t the  outer  scale,  S the  mean  strain  rate  (S  ~ njl),  £ the  rate  of  dissipation,  s the 
turbulent  strain  rate  and  using  6 1.5  f“  3k  , f ~ u f,  and  making  s/S  10, 

Tennekes  and  Lumley20  showed  that  isotropy  should  not  extend  down  much  further 
than  the  value  of  k where 


i 

I 


ki  350  . 

This  implies  that  the  maximum  wavelength  for  isotropy  would  not  be  much  greater 
than  about  2 percent  of  the  outer  lengt*n. 


17.  Pao,  Y.  H.  (1973)  Measurements  of  internal  waves  and  turbulence  in  two- 

dimensional  stratified  shear  flows,  Bound-Layer  Meteorol.  5:177. 

18.  Dutton,  J.  A.  (1970)  Effects  of  turbulence  on  aeronautical  sy  ms.  Progress 

in  Aeronautical  Sciences,  11,  edited  by  D.  Kuchemann  el  M,  Pergamon 
Press,  Oxford. 


19.  Ellison,  T.  H.  (1956)  Atmospheric  Turbulence,  Surveys  in  Mechanics 
(p.  425),  edited  by  T.  H.  Batchelor  and  R.M.  Davies,  Cambridge 
University  Press. 


20.  Tennekes,  H. , and  Lumley,  J.  L.  ( 1972)  A First  Course  in  Turbulence, 
Press,  Cambridge,  Massachusetts. 
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When  the  drain  of  kinetic  turbulent  energy  due  to  work  done  against  the  stable 
buoyancy  forces  is  taken  into  account,  Lumley"1  has  shown  that  the  inertial  range 
will  not  extend  to  lower  wavenumbers  than  k^  (inverse  "buoyancy  length")  given  bv 


(1) 


where 


N 


2 

1 


JL  £>0 

© dz 


(2) 


and  where  © is  the  average  potential  temperature,  g the  acceleration  of  gravity, 
and  Nj  is  the  buoyancy  frequency. 

In  the  HICAT  data  the  measured  thickness  of  turbulent  layers  ranged  from 

500  ft  to  7,  000  ft  with  most  of  them  under  3000  ft  or  about  1 km.  Taking  1 km 

as  the  outer  length,  the  first  criterion  above  implies  that  20  m would  be  the 

approximate  maximum  wavelength  for  isotropy.  The  bouyancy  length  criterion, 

according  to  Zimmerman  and  Loving,22  Monograph  I,  CIAP  gives  such  values 

as  15  m,  51  m,  and  40  m.  As  already  indicated,  the  data  showed  that  isotropy 

in  the  HICAT  measurements  extended  to  wavelengths  between  30  m and  45  m. 

Thus,  by  all  criteria  the  inertial  range  should  not  extend  much  beyond  about 

50  m or  so.  But,  as  already  mentioned,  the  -5/3  spectrum  (in  the  approximate 

sense)  extends  to  several  thousands  of  feet  if  not  tens  of  thousands  of  feet  in 

wavelength.  It  is  precisely  this  discrepancy  which  we  now  set  out  to  explain. 

23 

An  important  observation,  made  by  Gifford,  ~ is  that  one  should  expect  that 
spatial  aliasing  effects  will  significantly  extend  the  -5/3  range  of  wavenumbers  to 
smaller  values  in  the  one-dimensional  spectra  as  compared  to  the  three-dimen- 
sional spectrum.  In  the  example  studied  by  Gifford,  he  showed  that  the  -5  3 range 
of  the  one-dimensional  spectrum  extended  to  a value  of  k about  five  times  sr  aller 
than  the  smallest  -5/3  range  value  of  k of  the  three-dimensional  spectrum.  Thus, 
a three-dimensional  spectrum  with  inertial  range  extending  to  50  m could  appear 
out  to  250  m in  the  one -dimensional  spectruni  due  to  spatial  aliasing  effects.  But 
as  we  have  seen,  we  must  explain  -5  3 ranges  out  to  wavelengths  of  10,  000  to 
40,  000  ft  (say,  beyond  at  least  5 km).  This  could  only  happen  if  the  anomalous 

21.  Lumley,  J.  L.  (1904a)  The  spectrum  of  nearly  inertial  turbulence  in  a stably 

stratified  fluid,  J.Atm.Sci.  2 1 1 02. 

22.  Zimmerman,  S.  P. , and  Loving,  N.  (1975)  The  Xatural  Stratosphere  of  1 -74, 

CIAP  Monograph  1,  Final  Report,  DOT,  CIAP.,  DOT -TST-75-51 . 

23.  Gifford,  F,  (1959)  The  interpretation  of  meteorological  spectra  and  correla- 

tions, J.Meteorol.  10:344-340. 
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-5/3  range  in  the  three-dimensional  spectrum  extended  at  least  as  far  as  the 
"outer  length"  of  the  1 km  layer  thickness.  The  large  40,  000  ft  range  would  be 
compatible  with  an  outer  length  of  the  largest  layer  thickness  seen  in  HIC'AT.  In 
the  following,  we  therefore  set  out  to  explain  -5  3 ranges  all  the  way  to  outer 
length. 

An  important  consequence  of  the  anomalous  -5  range  is  that  it  can  lead  to 

24 

erroneous  estimates  of  f when  the  method  of  Stew  a r*  and  Grant  is  employed.  In 

other  words,  as  Pao1  has  written,  it  "...  should  be  considered  as  a caution  flag 

to  those  who  derive  turbulent  energy  dissipation  rate  i from  -.pectra  measured 

with  relatively  slow -response  instruments...".  Ibis  is  because  the  experimenter 

could  mistake  a -5/3  range  for  a true  inertial  subrange,  and  the  latter  might  begin 

* 

at  wavenumbers  higher  than  the  resolution  of  his  instruments. 

The  plan  of  this  paper  is  briefly  as  follows.  First  we  will  examine  the  equa- 
tions of  motion  as  well  as  the  spectral  equations  for  turbulent  motion  and  tempera- 
ture fluctuations.  The  physics  of  these  equations  and  their  individual  terms  will 
be  discussed,  and  a spectral  theory  of  K-H  turbulence  will  be  developed  which 
predicts  a -5/3  range  extending  from  the  Kolmogorov  microscale  to  the  integral 
or  outer  scale.  A simple  physical  picture  will  then  be  given  of  the  "dual  cascade" 
responsible  for  the  phenomenon. 


2.  'I'll  I HU  ii,»  UIOY'DI  Tl  Hill  l.l  M I.  SM)  U.MI’I  H\ll  hi: 
lilt  II  \ THINS  IN  STH  \ III  II  II  'III  Ml  IT  Oft 
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In  this  section,  we  shall  follow  the  treatment  given  by  Lin  et  al.  “ While 

more  suited  to  our  purposes  because  of  its  generality,  the  latter  treatment  is 

essentially  along  the  same  lines  as  that  given  in  such  standard  works  as  Hin?.e, 

27 

and  I.utnley  and  Panofskv.  “ First,  we  consider  the  Navier-Stokes  equation.  We 
assume  constant  viscosity,  incompressibility,  and  validity  of  the  Boussinesq 


A typical  relation  between  turbulent  diffusivity,  K and  the  dissipation  rate,  e, 
1/3  —4/3  • e 

is  K (const.)  e ' kj  w'here  N = buoyancy  frequency,  and 

^ -1/23/2  22 

kj  s (const.)  e ' N ' . (See  Zimmerman  and  Loving.  ) 

24.  Stewart,  R.W.,  and  Grant,  H.  L.  (1962)  Determination  of  the  rate  of  dissipa- 

tion of  turbulent  energy  near  the  sea  surface  in  the  presence  of  waves, 

J.  Geophys,  Res.  67:3177-3180. 

25.  Lin,  J.T.,  Panchev,  S. , and  Cermak,  J.  (1969a)  Turbulence  spectra  in  the 

buoyancy  subrange  of  thermally  stratified  shear  flows,  Project  Themis. 
Technical  Rept.  No.  1,  College  of  Engineering,  Colorado  State  t niversity. 

26.  Hinze,  J.O.  (1959)  Turbulence,  McGraw-Hill,  New  York. 

27.  Lumley,  J.L.,  and  Panofsky,  II.A.  (1964)  The  Structure  of  Atmospheric 

Turbulence.  Interscience  Publishers,  New  York. 
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; 


I II  III 


(7) 


r and  r.  x'  - x.,  local  homo- 

.1  _ .1  .1 

- • (dl '^/dx.j)  was  used,  and 

Uj  and  T W'ere  assumed  to  be  functions  of  only  the  vertical  coordinate,  x.?.  The 
reader  should  consult  the  references  for  further  details. 

We  now  physically  identify  terms,  and  since  they  all  involve  two  point  corre- 
lations, it  will  be  more  convenient  to  consider  the  zero  lag  case. 

I - rate  of  change  per  unit  mass  of  the  kinetic  energy  fluctuations  w ith 

respect  to  time. 

II  = rate  of  production  of  velocity  fluctuations  due  to  the  vertical  transport 

(u.?)  of  horizontal  fluctuations,  u^,  in  the  presence  of  vertical  mean 
shear,  dl^/dx.^. 

III  (as  will  be  discussed  below)  the  effects  of  the  distortion  and  rotation  of 
turbulent  motion  due  to  mean  shear. 

IV  - nonlinear  convective  effects  involving  vortex  stretching  of  smaller  eddies 

by  larger  eddies. 

Both  III  and  IV  are  divergences,  and  energy  is  neither  created  nor  destroyed  by 
them. 

V rate  of  work  done  against  stable  buoyancy  forces  which  results  in  the 
increase  of  the  potential  energy  of  the  temperature  (or  density)  fluctuation 
field. 

VI  = rate  of  turbulent  energy  dissipation  into  heat  by  viscosity. 

In  the  case  of  Eq.  7 we  have: 

I = rate  of  change  of  mean  square  temperature  fluctuations  with  respect  to 

time. 

II  = rate  of  production  of  temperature  fluctuations  due  to  vertical  transport 

(Ug)  of  temperature  fluctuations  (0)  in  the  presence  of  vertical  mean 
gradient  dT /dx^. 


d(0'0) 


dT 


9t  + <0'u3  + 0U3)  dx.5  + r3  ar^"  9x3 


- ( 0 1 G ) 


dU 


IV 


+ (00'u!  - 00'u!)  = 2v, 


32(0'0) 


9r . 
1 


w y <--T  gr  gr_ 

J J 3 


where  primes  refer  to  displaced  points  x'  = x + 
geneity  was  assumed,  the  approximation 
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III  = the  distortion  and  rotation  of  temperature  fluctuations  due  to  mean  shear 

which  is  analogous  to  term  III  in  the  equation  for  kinetic  energy  above. 

IV  = nonlinear  convective  effects  involving  vortex  stretching  as  in  the  kinetic 

energy  equation. 

V = rate  of  temperature  fluctuation  dissipation  by  molecular  diffusion  effects. 

The  spectral  equations  are  obtained  in  the  usual  manner  by  Fourier  transfor- 
mation of  the  two  point  correlations.  We  next  consider  the  three  dimensional 
scalar  wave  number  spectra  obtained  by  averaging  over  spherical  shells  in  k space. 
Equations  6 and  7 are  thus  transformed  into 


I 

9d(k.  t) 
8t 


II 

dU. 

i (k.  t)  ^ 
'uw  ’ dx0 


III 


9E^j 

T 9k„ 


SP.  AV. 


IV  V VI 

F(k,  t)  + ~ dwT( k,  t)  - 2v  k2  dfk,  t) 


and 


I 


II 


90TT(k,  t) 

at 


<WM) 


III 

aETi,i 

i ak3 


•SP.  AV. 


IV  V 

FTT(k,  t)  - 2isy  k2  dTT(k,  t) 


(8) 


(9) 


where  x*  s x3/2,  and  where  E.  . and  E,j.  . . are  the  diagonal  elements  of  the  non- 
averaged  three-dimensional  spectra. 

The  numbering  of  these  spherically  averaged  spectral  terms  corresponds  to 
the  previous  numbering  of  the  correlation  terms.  We  now  give  the  physical  inter- 
pretations. For  Eq.  (8)  we  have: 

I the  rate  of  change  with  respect  to  time  of  d(k,  t),  the  kinetic  energy 
spectral  density  at  k. 

II  = rate  of  " turbulent  production"  at  k.  <6  is  the  real  part  of  the  cross 

spectrum  between  the  streamwise  and  vertical  components  of  the  velocity 
fluctuations. 
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III  = transfer  in  k-space  of  kinetic  energy  due  to  the  distortion  effects  of 

mean  shear.  It  is  from  lower  to  higher  values  of  k in  the  anisotropic 

2 G 

region.  The  reader  may  consult  discussions  in  Hinze“  and  Lumley  and 
27 

Panofsky  for  further  information. 

IV  = transfer  in  k-space  of  kinetic  energy  due  to  the  inertial  cascade  from  low 

values  of  k where  production  is  important  to  high  values  of  k,  where 
molecular  dissipation  occurs.  It  represents  the  net  input  or  output  (that  is, 
the  "build-up")  of  energy  at  k and  when  F = 0,  the  cascade  is  steady  in 
k-space  since  F=  -3?(k)/3k. 

V rate  of  turbulent  buoyancy  " production,  " or,  as  in  the  present  stable 
case,  buoyant  "dissipation."  It  refers  to  kinetic  energy  loss  at  k by 
conversion  into  buoyant  potential  energy.  owT  is  the  real  part  of  the 
cross -spectrum  taken  between  fluctuations  of  temperature  and  vertical 
velocity. 

VI  = viscous  dissipation  rate  at  k considered  negligible  for  all  but  the  highest 

values  of  k near  the  Kolmogorov  microscale. 

The  terms  for  Eq.  (9)  are: 

I = rate  of  change  with  respect  to  time  of  bTT(k,t),the  temperature  fluctua- 

tion spectral  density  at  k. 

II  = rate  of  production  of  temperature  fluctuations  at  k.  Comparison  with 

term  V of  Eq.  (8)  shows  that  multiplication  of  the  present  term  by 
(g/T)/T-  (where  T'  s dT/dx*)  gives  term  V for  potential  energy.  In  fact, 
if  Eq.  (9)  is  multiplied  by  this  "conversion  factor,”  our  temperature 
fluctuation  equation  becomes  the  equation  for  the  potential  energy  of  the 
fluctuations. 

III  = transfer  in  k-space,  due  to  mean  shear.  Lin  et  al25  appear  to  be  the 

first  to  have  incorporated  this  term  in  the  equation  for  temperature 
fluctuations. 

IV  = the  inertial  cascade  term  for  the  temperature  fluctuations. 

V = molecular  heat  dissipation  rate  at  k which  is  assumed  negligible  for 

small  k. 

We  now  turn  to  a more  detailed  physical  description  of  the  terms  of  these 
equations  in  the  context  of  K-Il  turbulence. 
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3.  I Kinetic  hnrrsn  i ML  Produdimi  lerair- 

Insight  into  the  production  terms  II  and  V of  Eq.  ( f>)  in  the  context  of  K-II 

10  98  99 

turbulence,  is  provided  by  the  discussions  of  Ludlam  and  Businger”  ’ “ who 
derived  the  critical  Richardson's  numbers  on  the  basis  of  a simple  physical  model 
for  the  potential  and  kinetic  energy  budget.  In  the  following,  we  will  see  why  the 
production  of  turbulence  from  mean  shear  II  and  the  buoyancy  dissipation  (or 
negative  buoyant  "production"),  V,  are  approximately  equal  but  opposite  in  sign 
for  Kelvin-Helmholtz  turbulence.  Later  on,  this  crucial  observation  will  be  used 
to  explain  the  -5/3  spectrum. 

The  following  argument  leads  to  Richardson  number,  R;  0.25  as  a necessary 

condition  for  turbulent  breakdown  to  occur  when  there  is  a shear  across  a stable 
30 

layer  (cf.  also  Miles  ).  The  argument  is  based  on  the  fact  that  the  available 
kinetic  energy  from  the  shear  flow  must  exceed  the  potential  energy  of  buoyancy 
work  (involved  in  the  vertical  transport  of  parcels  of  fluid  across  the  layer  in  the 
mixing  process)  if  turbulence  is  to  take  place. 

Consider  the  idealized  case  of  a horizontal  layer  across  which  is  a linear, 
stable  temperature  gradient  and  a linear  shear  of  the  horizontal  velocity.  If 
parcels  of  fluid  of  unit  mass  at  the  top  and  bottom  of  the  layer  are  adiabatically 
interchanged,  one  must  do  work  against  the  buoyancy  force,  gA 0/0.  where  A0  is 
the  change  of  potential  temperature  across  the  layer  and  0 is  the  average  poten- 
tial  temperature.  A0  (d@/dx^)Ax3  where  Ax^  is  layer  thickness.  The  work 
done  to  move  one  parcel  across  the  layer  is 


Axl^  dAx!j 


g d0 

0 dx~ 
o 3 


(Ax3)2 

2 


(10) 


Previously  we  have  assumed  incompress 
the  more  general  case  treated  by  Businge. . 

28.  Businger,  J.A.  (1969a)  Note  on  the  critical  Richardson  number(s).  Quart. 

J.  Roy.  Meteorol.  Soc.  95:653 -654. 

29.  Businger,  J.A.  ( 1969b)  On  the  energy  supply  of  clear  air  turbulence,  in 

Clear  Air  Turbulence  and  its  Detection,  edited  by  Y.H.  Pao  and 
A.  Goldburg,  Plenum  Press. 

30.  Miles,  J.W.  ( 1963)  On  the  stability  of  heterogeneous  shear  flows.  Part  II, 

J.  Fluid  Mech.  16:209-227. 


ibility,  however  here  we  are  considering 
28‘ 
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The  potential  energy,  (PE),  generated  by  an  interchange  of  two  parcels  originally 
located  on  opposite  sides  of  the  layer  is  therefore 


PE  ( Ax.,)“ 

®o  dx3  3 


Next  consider  the  kinetic  energy,  <KE),  available  from  parcel  exchange  across 
the  layer,  due  to  the  shear.  Retting  L'j  be  the  mean  horizontal  velocity  and  Al  ^ 
be  the  difference  in  horizontal  velocity  across  t lit?  layer,  we  find  the  following 
available  KE  due  to  the  averaging  out  of  the  parcel  moments  after  interchange: 


KE  i pML1  + AU1)2-2(t1  + lAi:i)2] 


Equating  KE  PE  we  get: 

p oe/9x:j) 

R.  = (gradient  Richardson  number)  — — — — 

1 ®o  Ol  ^/9x.j) 


j 0.25  . 


We  now  relate  this  to  Eq.  ((!).  Retting  x'  x there,  we  see  that  term  II  is 
proportional  to  [u^u^  dEj/dx^J  while  term  V is  proportional  to  [u.^0  g T ] . The 
ratio  of  these  terms  is:  . 


R._  (flux  Richardson  number) 


T(UjUg)(dt  dx.^) 


Retting  K be  the  eddy  viscosity  and  the  eddy  conductivity  defined  by 


-u. u.,  = K dt  , / dx„ 
13  m 13 


tu.j  = dT 


we  obtain 


R.-SR.  when  K = I\u 
if  l m h 


In  the  case  where  intimate  mixing  can  occur,  KE  - PE  implies  R.  --  1 as 

29  1 

shown  by  Businger  with  arguments  similar  to  the  above.  This  in  turn  implies 

that  turbulence  will  begin  to  decay  when  R.  = 1 due  to  a lack  of  net  energy  input. 

In  the  case  of  K-H  turbulence,  the  value  of  R^  must  be  roughly  in  the  range  of 
0.  25  to  1.0  for  the  "active  life"  that  is,  the  energy  "fed"  period  of  the  turbulence 
layer.  It  should  be  now  clear  that  terms  II  and  V of  Eq.  (6)  physically  relate  to 
KE  and  PE  in  the  presence  of  the  velocity  shear  and  stable  potential  gradient. 

Next  we  consider  the  corresponding  terms  II  and  V of  the  spectral  equation, 
Eq.  (6),  and  we  see  that  they  are  the  spectral  form  of  the  KE  and  PE  at  a given  k. 
Their  ratio  is  the  spectral  flux  Richardson  number. 

Also,  note  that  the  "parcel  interchange"  picture  of  Businger  is  physically 
related  to  the  terms  0 u3  and  u^.  These  relate  to  the  PE  and  KE  of  the  turbu- 
lence due  to  Ug,  that  is,  vertical  transport.  I'nlike  Businger's  case,  incom- 
pressibility is  assumed  in  our  treatment;  thus,  potential  temperature  is  replaced 
by  ordinary  temperature. 

Another  point  is  that  the  production  of  turbulent  KE  can  also  be  regarded  from 

a point  of  view  that  does  not  involve  parcel  exchange  but  rather  an  energy  input 

from  mean  shear  to  turbulence  by  means  of  vortex  stretching.  As  Tennekes  and 
20 

Lumley  (P41)  showed,  a vertical  shear  tends  to  "feed"  vortices  by  this  process 
if  the  latter  are  aligned  at  45°  from  the  horizontal. 

i.2  I hr  Potential  1 fieri.  \ (I’f.f  I’rotltn  ttan  Irmi- 

As  mentioned  previously,  Eqs.  (7)  and  (9)  for  temperature  fluctuations  can  be 
regarded  as  the  equations  for  potential  energy.  For  example,  when  Eq.  (9)  is 
multiplied  by  (g/T)/T'  it  becomes  the  "PE  spectral  equation."  There  is  but  one 
"production  term,"  namely  II.  It  represents  the  rate  of  increase  of  the  squared 
temperature  fluctuations,  0 , in  the  presence  of  a mean  vertical  temperature 
gradient  though  the  action  of  vertical  transport  by  velocity  fluctuations.  For 
example,  an  increase  of  0 takes  place  when  a fluid  parcel  is  transported  down  the 
gradient. 

This  mechanism  of  production  is  entirely  analogous  to  the  momentum  case 
discussed  previously.  Also  notice  that  the  buoyancy  sink  term,  V of  the  KE  spec- 
trum equation,  Eq.  (8),  corresponds  to  the  buoyancy  source  term,  II  of  the  PE 
spectrum  equation,  Eq.  (9).  This  latter  point  explains  a certain  peculiarity 

pointed  out  in  discussions  of  the  so  called  "buovancy  subrange”  (see  for  example 
21  3 1 32 

Lumley,  and  Bolgiano,'  ’ ) namely  that  whenever  the  KE  spectrum  is  made  to 

31.  Bogliano,  R.  (1959)  Turbulent  spectra  in  stably  stratified  atmosphere, 

J.  Geophys.  Res.  154:2220. 

Bolgiano,  R.  (19(52)  Structure  of  tv.  ’•'uulenee  in  stratified  media,  J.Geophvs. 
Res.  67:3015-3023.  1 
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be  steeper  than  -5/3  by  a predominant  buoyancy  dissipation  effect,  the  Phi  spec- 
trum is  at  the  same  time  made  less  steep  than  -5/3.  These  slopes  bear  a recip- 
rocal relation  in  general,  and  the  above  description  enables  one  to  see  the 
physical  reason  for  this  effect. 
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As  we  have  seen,  the  nonlinear  convective  terms  of  the  equations  of  motion 
lead  to  the  divergence  of  triple  correlation  terms,  IV,  in  Eqs,  (6)  and  (7).  These 
give  rise  to  F and  FTT,  the  cascade  terms  in  Eqs.  (8)  and  (9).  In  an  inertial 
range,  these  terms  equal  zero,  implying  a steady  flow  down  the  scale  cascade  in 
k-space.  The  -5/3  slope  depends  on  a steady,  conservative  flow  (without  sources 
sinks,  "build-ups,'1  or  "decays"  along  the  way).  We  now  review  the  physics  of 
the  cascade  processes  in  the  PE  and  KE  spectra. 


1. 1 Inertial  I kb  I ( ascade 

Since  the  inertial  cascade  plays  a key  role  in  this  paper,  it  w'ill  be  useful  to 
examine  it  closely.  F’irst  wfe  will  consider  the  famous  treatment  by  Onsager33 
which,  to  date,  seems  to  be  the  most  informative  one  available.  Next  we  shall 
discuss  the  vortex  stretching  mechanism  responsible  for  it,  and  finally  the 
Heisenberg  concept  will  be  mentioned  since  it  gives  an  interesting  physical  picture 

and  will  be  used  later. 

33 

Onsager  proceeded  by  directly  Fourier  transforming  the  Navier-Stokes 
equations.  The  nonlinear  convective  term  depends  upon  the  difference  of  wave- 
numbers,  (t  k1  ± k).  Assuming  that  k,  as  well  as  differences  of  k,  are  of  the 
order  of  (1/L)  where  "L"  is  the  size  of  the  largest  eddies,  he  found  that  k'  would 
be  at  most  2/L,  where  k is  regarded  as  the  driving  mode  and  k'  the  mode  receiving 
the  energy.  Applying  this  reasoning  to  subsequent  cascade  steps  results  in  a geo- 
metric progression.  For  this  reason,  the  wavenumber  doubles  at  each  step; 
hence,  Ak  between  steps  ~ k.  The  amount  of  energy  transferred  at  each  step  is 
therefore  of  the  order  kb(k).  Let  Q(k)  represent  the  flux  of  energy  past  wavenumber 
k in  k-space.  Then,  in  the  continuous  limit, 


Q(k)  = 


k 6(k) 
r(k) 


(lfi) 


33.  Onsager,  L.  ( 1949)  Statistical  hydrodynamics,  Nuovo  Cimento  f>:(Ser.  9 
Supp.  2):279-287.  ’ 
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r 


r 


U 


where  r(k)  is  a characteristic  time  for  a step  at  wavenumber  k.  On  apparently 
dimensional  grounds,  Onsager  chose 


T 1/V 


k 0<k>  . 


34 

Alternatively,  Corrsin'  deduced  this  result  physically  from 


(17) 


kinetic  energy  per  step 
7 energy  transfer  rate 


(18) 


He  obtained  a denominator  by  analogy  with  u^u^  (3  U./3x^): 

3/2 


(spectral  velocity)'  (k <t>) 

(spectral  length)  (1/k) 


V 


i 5.3 

k o 


(19) 


Hence 


r(k) 


k<£ 


y k5b3  ^ k3p' 


QED. 


Using  this  in  Eq.  (16) 
Q(k)  * b(k)3/<2  k5f2 


(20) 


(21) 


Assuming  that  the  cascade  in  k space  is  divergenceless  and  conservative. 


3Q 

75TT 

Setting: 


= -F  = 0 


(22) 


Q(k)  e 
<J(k)  * e 


. .2/3  k — 5/3 


(23) 

(24) 


This  of  course  agrees  with  the  usual  result  obtained  by  a similitude  argument. 

The  doubling  at  each  step  accelerates  the  cascade  and  therefore  causes  the 
amount  of  energy  at  k to  decrease  with  increasing  k in  a way  which  is  analogous  to 
the  manner  that  a convectively  accelerated  stream  of  incompressible  fluid  will 
become  more  narrow  as  it  speeds  up. 


34.  Corrsin,  S.  (1958)  Local  Isotropy  in  Turbulent  Shear  Flow,  NACA  KM  58-B11. 
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Physically,  the  cascade  comes  about  by  the  vortex  stretching  action  of  the 
larger  eddies  on  the  smaller  ones.  Mathematically,  this  is  represented  by  term 
IV  in  Eq.  (C).  It  can  be  shown,  however,  that 


3 I 1 \ 

ST—  — U-U-U.  U.U.S.. 

9x.  I 2 113/  1 .1  13 


(2  5) 


where 


, , 3u.  3u. 

1 I L 3 


-ij  2 3 xj 


(26) 


But  (-u.u.s.  ) is  the  amount  of  energy  per  unit  mass  and  per  unit  time  gained  b\  a 
1 .1  l3 

disturbance  with  velocity  components  u.u.  in  strain  rate  s..  (see  Tennekes  and 
OH  1 J •' 

Lumley  p.  257).  If,  for  example,  there  is  compression  in  the  y direction  and 

stretch  in  the  x direction  (Figure  2)  the  energy  exchange  rate  is  equal  to 


Q = s(u;  - u“) 


(27) 


Since  ^ is  increased  and  ^ decreased,  u!J  increases  while  u^  decreases.  Hence, 
(u^  - u”)  starts  at  zero  and  becomes  positive.  The  not  effect  is  to  increase  the 
energies  of  eddies  bv  this  deformation.  The  cascade  is  therefore  visualized  as 
the  transfer  of  energy  by  the  distortion  effects  of  large  eddies  upon  the  smaller 
eddies.  The  interested  reader  should  consult  the  excellent  discussion  of  Tennekes 
and  Burnley‘S  for  further  quantitative  information  about  this  process.  They  show 
for  example  that  this  leads  to  a derivation  of  the  -5  3 slope.  Since  v ~ ko 
dimensionally. 


'■  1 
I 

, i 


Q(k)  ~ kq(k)  s(k)  , 

and,  also  dimensionally 
'1  1/2 

s(k)  ~ <k  $)  . 

Equations  (23),  (28),  and  (29)  thus  lead  directly  to 

p(k)  ~ e2/3  k‘5/3  . [see  Eq.  (24)] 

The  computationally  useful  concept  of  Heisenberg  is  that  the  energy  transfer 
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across  k can  be  regarded  as  a turbulent  viscosity  effect  (see  Ilinze  ): 


(28) 


(29) 


/ 


F dk 


’jf 


2 r(k)  / k"  b(k)  dk 


(30) 


where  r(k)  is  the  effective  turbulent  viscosity  defined  (from  dimensions) 


where  a ~ 1. 

Various  generalizations  of  this  approach  have  been  advanced  and  we  shall 

3 5 

later  make  use  of  a very  general  one  due  to  Panchev. 

1.2  I rmiHTHturi'  I Im  liialton  PI  I ascailr 

The  mechanism  of  the  previous  inertial  cascade  represents  a small  version 
of  the  ' turbulent  production  process'  due  to  average  shear.  In  a similar  manner, 
the  temperature  cascade  is  a small  version  of  the  production  process  due  to  mean 
temperature  gradient. 

As  above,  the  cascade  term  IV  in  Eq.  (7)  is  the  divergence  of  a triple  corre- 
lation. Considering  IV  as  a one  point  correlation,  we  follow  the  treatment  of 
Tennekes  and  I.umlev.  Briefly,  IV  can  be  represented  as  a miniaturized  form  of 


35.  Panchev,  S.  (1971)  Random  Functions  and  Turbulence,  Pergamon  Press, 
New  York. 
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( Oil 


m 


(temperature  \rariance  production  rate) 


[32) 


Ue  estimate  nu.,  by  (9T  9x^1  and  where  n and  ( are  representative  velocity  and 
length. 


^ k,^  o(k.j) 


f ~ k.;1  . 


(33) 


(34) 


The  miniaturized  (9T /9x.^),  that  is,  the  gradient  of  the  "larger  eddy"  acting  in  the 
"smaller  eddv"  is 


(gradient  of  temp) 


0 rj-  ( k ) 


(35) 


k.,  > k;.  I.et  Q be  the  flux  past  k in  the  temperature  cascade.  Wo  than  have 

Q0  = k2l  V k2  0<k2)l  kl  °TT(k)  • 

Assuming  kj  ~ k., 

c k’’  “ o-j.rj-(k)  o(k)  ' . 


( :i  fi) 


(37) 


This  derivation  shows  the  important  point  that  the  scaler  temperature  cascade 
is  due  directly  to  the  action  of  the  velocity  field  on  the  temperature  field.  This 
important  point  will  be  used  below. 
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The  only  way  to  explain  an  anomalous  -5  i spectrum  is  to  assume  that  one  has 
a situation  which  is  equivalent  to  the  conservative  steady  cascade.  In  the  present 
case  we  argue  that  in  K-H  turbulence,  terms  II  and  V in  Eq.  ( 8)  essentially  cancel 
each  other  over  all  wavelengths  from  microscale  to  outer  scale.  This  can  be  true 
only  if  the  two  terms  are  approximately  equal  and  opposite  and  in  addition  have  the 
same  dependence  on  k.  To  prove  the  latter,  it  is  necessary  to  show  that  the  tem- 
perature field  is  determined  by  the  velocity  field.  This,  however,  has  already- 
been  done  by  I,umley“'  and  the  necessary  assumptions  are:  (1)  both  Reynolds 
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number  and  Peclet  number  must  be  high  enough  that  molecular  effects  can  be 
ignored.  (2)  OT/3t)  = 0;  and  (3)  {sr2  a2T/3x2  Le}/OT79x.)«  1 where  L^.  is  the 
Eulerian  space  integral  scale,  L„  ~u'T.  (Corrsin2^)  where  u'  is  the  r.m.  s. 
velocity  and  Tj  the  Lagrangian  time  scale.  Strictly  speaking,  the  mean  tempera- 
ture is  changed  in  time  in  K-H  turbulence.  However,  we  make  the  approximation, 
for  an  interval  of  interest,  that  we  have  stationarity. 

Briefly,  Lumley's  proof  begins  with 


M+  u (&-  + 
at  i\  axt  ax.  ) 


= o 


(38) 


which  is  Eq.  (5)  with /jq-  =0,  T = T + 6,u^  = Ub  + m and  3T/3t  - 0.  Equation  (38) 
is  in  Eulerian  form.  Burnley  solved  (38)  under  his  above  stated  assumptions  to 
arrive  at 


9(x,  t) 


(39) 


where  a(x,  t)  is  the  position  at  t ■-  0 of  that  fluid  particle  which  will  arrive  at  x at 
time  t.  In  our  problem,  where  we  assume  a vertical  temperature  gradient  con- 
stant in  space  and  time,  we  have  (OT/ax^)  = T') 


i9(x,  t)  T’(z  - z ) 
o 


(40) 


Physically,  Eq.  (40)  is  easy  to  understand.  The  temperature  fluctuation  is 
all  due  to  the  vertical  displacement  of  parcels  from  their  initial  position  to  new 
positions  where,  by  virtue  of  the  gradient,  they  are  at  temperatures  different  from 
the  mean  (compare  Section  4.2  above). 

The  conclusion  of  Lumley's  proof  is  that,  from  Eq.  (39)  we  can  infer 


0(x,  t)  u.  (x,  t) 


(x.  - a^  u..  (x1,  t) 


s (-aT/aXjMA..)  . 


(41) 


Since  A.  depends  only  on  the  velocity  field,  the  proof  is  completed. 

9 1 

Using  the  fact  that  the  velocity  field  controls  the  temperature  field  Lumley“ 
went  on  to  show,  by  means  of  a similitude  argument  that 


*wT 


-C 


a T 

1 3x.^ 


1 / 3 . -7/3 
e k ' 


(42) 


36.  Corrsin,  S.  (1963)  Estimates  of  the  relations  between  Eulerian  and  Lagrangian 
scales  in  large  Reynolds  number  turbulence,  J.Atm.Sci.  20:113-119. 
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where  C ^ is  a constant  of  order  one.  This  result  is  an  application  of  the 
Kovasney37  formalism  which  says  the  spectral  forms  are  determined  by  the  local 
spectral  flux. 

It  should  be  mentioned  that  l.umley's  proof  of  Eq.  (4  2)  rested  on  one  more 

3 8 

assumption,  namely  "local  inertialitv"  (cf.  l.umley  ): 

(e/k)/OE/ak)»  1 . (43) 


Physically,  this  assumption  means  that  the  I energy  in  or  out  of  the 

cascade  is  but  a small  fraction  of  the  throughput.  I hi-  assumption  is  consistent 
with  the  conclusion  which  follows  from  it  I el  .v. 

In  a similar  way,  one  can  show 


6 


uw 


(44) 


where,  hopefully,  C ^ = C.,. 

in  Eq.  (i!)  let  us  assume  that  at  the  value  of  k interest,  we  can  ignore 
molecular  dissipation;  that  except  for  the  effects  of  produr'ion  terms,  the  cascade 
is  steady  and  conservative;  and  that  q ~ 0.  In  tnis  case,  using  I - (-de  '3k), 


9f 
0 k 


dll, 

f>  — — + -M-  A 
’’uw  dXg  — °wT 


e1/3  k-7 


C(1  - Hi) 


(45) 


where  9«/9k  is  the  change  of  throughput  with  respect  to  k,  and  Hi  is  the  gradient 
Richardson  number.  Recalling  Hi  s [(g  f)-dT  dx..,|  (dl'j  dx^P  we  see 


0.  25  < Ri  < 1 


(40) 


for  K-II  turbulence  in  view  of  the  previous  discussion.  This  implies  that  (9e  (9k) 

will  remain  small  except  for  large  wavelengths.  In  order  that  Oe  9k)  e k,  we 

have  as  a lower  limit  on  k,  k , 

c’ 


It  is  assumed  that  Ri  Rf. 

57.  Kovasznay,  E.  S.  ( 1P4H)  Spectrum  of  locally  isotropic  turbulence,  J.  Aeron. 
Sci.  15:741-753. 

38.  l.umley,  J.  E.  (lfl(>7a)  Theoretical  aspects  of  research  on  turbulence  in 

stratified  flows,  in  Atmospheric  Turbulence  and  Radio  Wave  Propagation. 
edited  by  A.  M.  Yaglom,  and  V.  I.  Tatarsky,  Publishing  House  "Nauka,  " 
Moscow. 
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kt:7/3(U\)(l  - Hi)  e — 


where 


(Uj)  = dUj/dXg 


The  following  estimates  will  be  employed: 


(U’j)  ~ U/L 


where  L is  the  outer  length.  Also 


E ~ 77  • Cj  (50) 

where  u is  the  rms  value  of  the  velocity  fluctuations  and  C(  is  a constant  which 
takes  the  buoyancy  into  account. 

12 

According  to  Woods  and  \\  iley,  ” the  rate  at  which  turbulence  loses  energy 
to  viscous  dissipation  is  much  greater  (up  to  a factor  of  10)  than  the  rate  at  which 
the  turbulence  does  work  against  gravity  by  changing  the  density  (temperature) 
profile.  This  implies  is  near  unity. 

Defining  C2  as 

u c 2 1 V 1 1 (51) 

and  Eqs.  (49),  (50)  and  (51)  in  (47),  we  obtain 


.4/3  <1-Ri) 

C2/3C2 

1 2 


Letting  Ac  be  the  maximum  wavelength  for  inertiality,  that  is,  k,  = 2y/ Ac,  we  have 

1 /o  3 /2  t /4 

Ac  (2tLC|' -C2'2)/(l  - Ri)3  4 . (53) 


From  this  relation  it  follows  that,  as  long  as  Cj  and  C„  are  not  zero  (that  is, 
assuming  there  is  turbulence  and  thus  u * 0)  there  always  exists  a value  of  Ri  less 
than  unity  such  that  Ac  L.  Thus,  the  -5/3  range  can  extend  to  the  outer  length 
for  Ri  sufficiently  close  to  unity.  The  effect  described  by  Gifford23  will  make  the 
-5/3  range  extent  to  about  5 Ac  in  the  one-dimensional  spectrum. 
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In  addition  to  all  the  above  mentioned  explicit  assumptions,  it  was  also 
implicitly  assumed  that  the  dynamic  effects  of  the  inertial  cascade  are  not  small 
in  their  effects  in  comparison  to  mean  shear  effects,  in  other  words  we  assume 
that  the  mean  strain  rate  is  not  large  compared  to  the  eddy  strain  rates.  In  our 
judgment,  all  the  above  assumptions  are  appropriate  for  the  approximate  descrip- 
tion of  K-H  turbulence.  This  weak  interaction  assumption  is  justified  on  the  basis 
that,  unlike  boundary  layer  or  wall  turbulence,  the  mean  structure  is  too  weak  to 
permanently  "feed"  the  turbulence  directly.  In  actual  fact,  the  layer  is  fed,  to  a 
large  degree,  at  the  turbulent -laminar  interface,  bv  laver  spreading 
(c.f.  Woods8'7). 
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We  now  put  the  spectral  equations  into  a form  more  amenable  to  their  resolu- 
tion and  visualization.  First  we  integrate  Eqs.  (8)  and  (9)  from  0 to  k and  then 
integrate  them  from  0 to  « and  subtract  the  former  from  the  latter,  t sing 

00 

£ J 2i'k^tf  dk  (KE  dissipation  rate)  , ... 


/2l'Tk2,tTT  dk  (temperature  fluctuation  (PE)  dissipation  rate) 


(Winze,  * p.  179),  we  obtain 

1 11  HI  IV  V VI 

* k k - * 1 
/ d dk  -£  + 2r/k%dk  -f  F dk  *4-  f 6 dk  . p f 0 

k ^ * q ' J 1 J WU 


(ki 


3 /SP.AV. 


for  the  KE  spectrum,  and 


f dTT  dk  <-\)  + 2 vT  j k2  dTT  dk  - J FTT  dk  - T'  J </>wT  dk 


- r I,  iV 

1 J ;’k3 


SP..AV. 


for  the  PE  spectral  equation,  where  we  have  used 


-J* 


SP.  AV. 


, / 9E-. 

'/  (ki  akT" 


Referring  back  to  the  identification  of  terms  for  Eqs.  (8)  and  (('),  the  physical 
meaning  of  most  of  the  above  is  straightforward.  As  for  the  "k-transfer'  terms 
IV  and  VII  in  Eq.  (56)  and  IV  and  VI  in  Eq.  (57),  they  refer  to  the  energy  (PE  or 
KE)  transferred  out  from  the  range  0 to  k into  the  range  k and  higher. 

In  the  following,  we  shall  consider  Eq.  (57)  as  having  been  multiplied  by 
(g,/T)(T')~*  so  that  we  can  regard  it  literally  as  the  PE  spectrum.  Equations  (56) 
and  (57)  thus  describe  two  energ\  cascades  that  are  coupled  through  the  (except 
for  sign)  common  term  V.  Figure  1 giver,  a schematic  diagram  of  the  dual  cascade 
represented  by  these  equations.  We  imagine  the  cascades  as  being  analogous  to 
pipes,  large  at  the  large  scale  (or  small  k)  end,  (symbolizing  the  large  amount  of 
energy  stored  there),  and  small  at  the  high  k end.  For  the  most  part,  the  "narrow  - 
ing of  the  pipes"  is  due  to  the  acceleration  of  the  cascade  discussed  previously. 

But,  at  the  high  k end  we  imagine  the  pipes  as  becoming  increasingly  porous  to 
symbolize  the  dissipation  from  the  molecular  effects,  \'  and  t.  The  process  of 
"going  down  the  pipes"  is  the  cascade. 


Woods  and  Wiley*2  point  out  the  e is  significantly  larger  than  X.  This  is  clear 
from  the  fact  that  Ri  is  presumably  below  0.  25  at  the  start  of  the  turbulence. 
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Figure  3.  Dual  Cascade  Model 


The  source,  it  we  ignore  the  <*i  term,  is  the  turbulent  production  at  the  large 
end  of  the  KE  cascade.  We  imagine  the  energy  as  entering  through  the  wall  of  the 
pipe.  The  PE  cascade  is  fed  by  draining  from  the  KE  cascade  through  the  "porous 
pipe  walls."  This  can  be  considered  as  a generalization  based  upon  Lumley's^1 
"porous  pipe"  concept.  Note  that  between  the  "fed"  scale  and  dissipation  scale 
there  is  a conservative  region  for  both  pipes,  and  this  of  course  represents  the 
inertial  range. 

The  most  important  point  to  note  about  the  above  picture  is  that  it  suggests 
the  physical  explanation  of  the  anomalous  -5/3  range.  If  we  imagine  that,  in  the 
KE  cascade,  the  input  due  to  production  and  output  due  to  "buoyancy  drain"  are 
approximately  equal  to  each  other  for  all  values  of  k where  they  are  important, 
then  their  net  effects  will  approximately  compensate  or  cancel,  and  the  "shape  of 
the  KE  pipe,"  that  is,  the  slope  of  the  turbulence  spectrum,  must  approximate  the 
-5/3  shape.  This  is  because  the  "flow"  would  then  only  depend  on  the  cascade 
term,  F,  which  accelerates,  as  Onsager  demonstrated,  in  just  such  a manner  as 
to  give  -5/3  as  slope. 

The  above  explanation  leaves  out  a number  of  things.  In  the  first  place,  tiie 
direct  spectral  transfer  due  to  mean  shear  has  been  omitted.  More  importantly, 
if  production  and  buoyancy  terms  exactly  cancelled,  the  KE  spectrum  would  have 
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no  source  at  all.  In  K-H  turbulence,  however,  if  we  assume  Ri  0.  25  for  insta- 
bility to  occur,  then  we  are  guaranteed  that  during  the  initial  phase  of  the  turbu- 
lence, production  will  exceed  "drainage”  so  that  there  will  be  a net  feeding  f the 
"KE  pipe,"  and  relation  (43)  assures  us  that  "local  inertiality”  will  exist  and  hern  >- 
to  a good  approximation  the  -5/3  shape  will  prevail.  As  for  mean  shear  effect - 
on  the  cascade,  we  have  already  assumed  that  the  interaction  between  mean 
and  turbulence  vorticity  is  "weak"  in  K-II  turbulence,  and,  even  including  , he  . r 
transfer"  effects,  it  will  be  shown  below  that  we  would  still  expect  - ■ .r  t he 

slope  in  the  "anisotropic"  region. 

Once  Ri  = 1 and  the  decay  phase  of  the  turbulence  begins,  the  "source"  f »h«- 

turbulence  will  be  the  slowly  decaying  large  eddies  (d)  and,  since  both  shear  and 

temperature  gradient  within  the  layer  will  have  been  destroyed  by  the  mixing 

process,  we  would  expect  a -5/3  slope  on  the  basis  that  local  inertiality  would  hold 

even  more  accurately  than  during  the  initial  phase.  In  this  second  phase,  both  K1 

and  PE  cascade  should  have  a -5/3  slope  which  is  in  contrast  to  the  initial  phase 

where  the  PE  slope,  in  the  anisotropic  region,  would  be  significantly  flattened  (if 

not  in  fact  reversed  in  sign)  due  to  the  effects  of  " production"  without  a compensat- 
*!e 

ing  drain.  On  the  other  hand,  at  the  very  start  of  the  turbulence  regime,  we 
might  expect  a k 1 dependence  in  the  KE  spectrum  due  to  "strong  interaction" 

(cf.  Gisina^). 

<1.2  The  \nomalous  -.'>  Range 

We  now  proceed  in  the  spirit  of  Section  5 to  derive  the  KE  spectrum  in  the 
2 1 

manner  used  by  Lumley  in  his  derivation  of  the  buoyancy  range  equations.  We 
start  with  Eq.  (45)  and  integrate  to  obtain 

,V3,,2/3  (1.£,iT1,2(1.R11,;S/3k-4/3) 

where  the  constant  of  integration  has  been  chosen  to  make  s when  k — on 

(t  is  the  viscous  dissipation).  Now  the  condition  given  by  Eq.  (43)  assures  us 
that  the  spectrum  will  be  determined  locally  by  the  spectral  flux  and  that  the 

A log-log  plot  of  the  temperature  spectrum  derived  from  HICAT  showed 
(Vinnichenko  and  Dutton^0)  a -5/3  power  law  on  the  average.  The  implication  of 
the  above  reasoning  is  that  most  of  the  data  were  taken  in  the  decay  phase. 

39.  Gisina,  F.A.  (1966)  The  effect  of  mean  velocity  and  thermal  gradients  on  the 

spectral  characteristics  of  turbulence,  Izvest.Atm.  and  Oceanic  Phvs. 
2:804-813. 

40.  Vinnichenoko,  N.  K. , and  Dutton,  J.A.  (1969)  Empirical  studies  of  atmos- 

pheric structure  and  spectra  in  the  free  atmosphere.  Radio  Sci.  4:1115- 
1126.  
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spectrum  will  remain  approximately  isotropic  since  anisotropic  inputs  or  outputs 
are  a small  fraction  of  the  "throughput"  e.  In  other  words,  we  have  the  condition 
where  the  cascade  is  "locally  inertial  in  wavenumber  space"  (a  phase  attributed  to 
Corrsin  by  Lumley).  This  allows  the  use  of  all  the  customary  forms  familiar  in 


the  inertial  subrange,  but  with  "variable  e"  as  Lumley 
Thus  we  find  that 
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has  put  it. 


<f>(k) 


-5/3  £2/3  [1.C(C,)2(1 


.-2/3 


-4 

Ri)e  ‘•/'J  k /3 
o 


(61) 


From  this  we  immediately  arrive  at  the  conclusion  that:  (1)  for  high  enough  k we 
have  the  ordinary  inertial  range  relationships,  and  (2)  under  the  restriction  of 
Eqs.  (43),  the  -5/3  relation  holds  for  all  scales  up  to  the  outer  length  which  is 
also  indicated  by  Eq.  (47)  whenever  we  are  concerned  with  K-H  turbulence  at  a 


stage  where  Ri  is  sufficiently  close  to  unity.  It  should  be  mentioned  in  passing 
4 1 

that  Lumley  also  explained  an  anomalous  -5/3  spectrum,  in  a different  context. 


by  means  of  this  constant  flux  mechanism. 

Another  approach  to  the  anomolous  -5/3  slope  is  the  following.  Vh.j  time  we 
use  Eq.  (56)  to  derive  the  -5/3  slope  of  the  KE  spectrum  of  K-H  turbulence  under 
the  assumptions  of  weak  interaction  between  turbulence  and  mean  profile,  and 
stationarity  (of  both  mean  and  turbulent  quantities).  We  shall  make  use  of  the 
Heisenberg  concept  mentioned  previously.  In  addition,  we  shall  make  use  of 
Tchen's  idea  of  using  the  Boussinesq's  concept  to  represent  terms  VI  and  VII 


together  as  a "dissipative  effect"  (see  Tchen^3  and  Hinze3^).  We  shall  also  treat 


term  V in  a similar  way  (cf  Gisina39). 
Thus,  for  terms  VI  and  VII  we  have: 


U' 
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d E . . % | _ « 

«wu-  kl  Tkr)ldk  -‘m(k><r'> 
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where 


rm(k)  a 


*•  / V 


<Mk) 


kJ 


dk  o'a  K (k) 


(63) 


41.  Lumley,  J.  L.  (1967b)  The  inertial  subrange  in  nonequilibrium  turbulence, 
in  Atmospheric  Turbulence  and  Radio  Wave  Propagation,  edited  bv  A.  M. 
Yaglom  and  V.  I.  Tatarsky,  Nauka,  Moscow. 


42.  Tchen,  C.M.  (1953)  On  the  spectrum  of  energy  in  turbulent  shear  flow, 
J.  or  Res.  of  NBS  50:51  -62. 
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and  for  term  V 


,Hik))T' 


r^j(k)  a"  a J +1  dk  »"o'K(k)  . 

k * k 


We  can  thus  write  Eq.  (56)  as 


e 2(i/  + eK(k))  J k2  6 d k + aK(k) 


a = a'o  (V  )2  - Tp-e"  <> T 


To  solve  (66),  one  uses  the  technique  of  Agostini  and  Hass  and 
Chandresekhar^  namely,  that  of  defining 


r(k)  = 


/ *2 


d(k)  dk 


1 dr 


In  the  resulting  differential  equation,  which  can  be  solved  by  direct  quadrature, 
the  boundary  condition  is  r ■ e/2r  when  k » « , giving  for  constant  of  integration 


co*ST  it  (-rhz 7)  • 


43.  Agostini,  L. , and  Bass,  J.  (1950)  Publ.  Sci.  et  Tech.  Ministere  Air  237. 

44.  Chandrasekhar,  S.  (1949)  On  Heisenberg's  elementary  theory  of  turbulence, 

Proc.  Roy. Soc.  London. 200A:20. 
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Thus,  as  Gisina^y  has  shown  in  a different  context, 


*<k)  = 


(e  + av) 


(3 am  (e  + av) 


For  values  of  k of  interest,  this  gives  a k ' dependence.  The  most  impor- 
tant part  of  this  result  to  notice,  other  than  the  -5/3  slope,  is  the  fact  that  the 
"constant"  no  longer  simply  depends  on  e,  as  it  would  in  a classical  inertial  range. 
Rather,  it  also  contains  (av)  which  involves  the  production  and  bouyant  dissipation 
terms.  As  mentioned  in  the  introduction,  this  could  lead  to  an  erroneous  deter- 
mination of  e if  the  investigator  makes  the  mistake  of  perceiving  an  anomolous  -5/3 
range  as  an  inertial  range. 

It  should  not  be  overlooked  that  Eqs.  (61)  and  (71)  are  not  consistent.  Only 

the  -5/3  exponent  is  consistent.  This  lack  of  consistency  is  excusable,  perhaps, 

on  the  grounds  that  the  Heisenberg  approach  usually  does  give  physically  incorrect 

2 fi 

results  outside  of  the  wavenumber  range  of  primary  interest  (cf  Hinze’s 
discussion  of  the  Heisenberg  result  for  the  dissipation  range).  * In  view  of  the 
necessity  for  an  inertial  range  between  the  anisotropic  production  range  and  the 
dissipation  range,  we  should  choose  the  relation  (61)  as  being  more  correct  at 
large  \ since  it  includes  the  inertial  range  as  a high  wavenumber  limit  in  contrast 
to  (71)  w'hich  does  not.  From  (61),  the  remark  above  about  determination  of  c 
still  holds. 


7.  -IM/CTHM.  \<>M  INI. \UI  I II  - 

In  both  HICAT  data  (Crooks  et  al4>)  and  in  the  tropospheric  data  of  Reiter  and 
Burns4*’  it  has  been  noted  that  there  are  "humps"  and  "dips",  sometimes 
called  "gaps"  in  the  KE  spectra  Various  attempts  have  been  made  to 
explain  these  (Cf  Reiter  and  Burns4'’)  One  such  explanation,  Pao4  ' . was 
based  on  a numerical  simulation  in  which  it  was  found  that,  with  high  degrees 


*The  explanation  for  the  somewhat  incorrect  nature  of  the  Heisenberg  formalism 
is  that  it  reflects  distant  wavenumber  interactions,  whereas  close  wavenumber 
interactions  are  in  reality  much  more  important. 

45.  Crooks,  W.M.,  Hoblit,  F.  M, , and  Prophet,  D.T.,  et  al  ( 1967)  Project 

HICAT:  An  Investigation  of  High  Altitude  Clear  Air  Turbulence,  Tech. 

Rept.  AFFDI ,-TR -67 -123.  AD824  865,  AD824  904,  AD825  369  (1967). 

See  also  AD846  086,  AD847  497  (1968). 

46.  Reiter,  E.R.  and  Burns,  A.  (1966)  The  structure  of  clear-air  turbulence 

derived  from  "TOPCAT"  aircraft  measurements,  J.  Atm.  Sci.  23:206-212. 

47.  Pao,  V.  H.  (1968)  Turbulence  Velocity  and  Scalar  Spectra  in  Stably  Stratified 

Fluids,  Boeing  Scientific  Research  Laboratories,  D 1-82-0680,  AD677-585. 
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of  stability,  the  nonlinearities  would  occur.  I'nfortunately  in  these  simulations  the 
production  term  due  to  shear  was  omitted. 

Another  possible  explanation  is  the  following.  In  Eq.  < 5 G),  let  us  ignore 
terms  I,  III,  and  VII,  and  in'Eq.  (57)  omit  I,  III,  and  VI.  We  then  have 


II 
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(72) 
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(73) 

25  48 

These  equations  have  been  treated  by  Lin  et  al“‘  ’ and  by  taking  the  limiting 
form  of  the  generalized  eddy-viscosity  approximation  due  to  Panchev'^  they 
arrived  at  the  following  nondimensionalized  representation 


IV 

C5/2  <l>  (X)3 /2  ! 


1 ( M 

\ 2 / 

r,  x v 

l 

1 '2 


C4/2 


<MX)  <t>TT(X) 


1 

j,  , (-i)  + (-^iL)  (l/2)+(C2,'2)J 

+ j ' r | x “ x 4>  (X) 

V 

v5/2  .1/2 
X ' <t> 


4>tt(X) 


1^4  C 2 

- -r-  + 1/2  4 ' - 

!rr | x 2 2 <t>(x)  4>tt(x) 


(74) 


(7.5) 


* 

where  Eq.  (74)  corresponds  to  (72),  and  (75)  corresponds  to  (73).  The  numbering 
of  terms  is  retained.  Also,  <t>  and  4'^.^.  are  the  nondimensionalized  spectra,  and 
X is  the  nondimensionalized  wavenumber.  1 a parameter  characterizing  tur- 
bulence production  due  to  mean  shear,  F.  a parameter  characterizing  the 
bouyancy  dissipation  which  depends  on  mean  temperature  gradient,  and  I'j,  = a 
parameter  characterizing  the  production  of  temperature  fluctuations  which  depends 
temperature  gradient  Parameters  CL  and  ('4,  both  of  which  can  range  from  0 to  1, 
represent  a generalization  of  Tchen’s  concept  of  strong  and  weak  interaction  between 
turbulence  and  mean  shear  and  temperature  gradient  respectively.  "0"  means 
weak  interaction  and  "1"  means  strong  interaction  in  the  sense  that  the  vorticitv 

Compare  term  IV  of  Eq  (75)  with  Eq  (37) 

48.  Lin,  J.T.,  Panchev,  S. , and  Cermak,  J.E.  (19fi9b)  A modified  hypothesis 
on  turbulence  spectra  in  the  bouyancy  subrange  of  stably  stratified  shear 
flow,  Kadio  Sci.  4:1333-1337. 
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squared  term  (in  the  Heisenberg  expression  for  ' dissipation'  ) is  the  square  mean 

gradient  or  shear  in  the  "weak"  case,  and  the  square  vortieity1'  or  equivalent  form 

26 

for  temperature  fluctuations  in  the  range  0 to  k in  the  "strong1'  case  (cf  Hinze 
pp.  264,  265).  Finally,  note  that  in  Eqs.  (74)  and  (75),  the  spectral  transfer  terms 
are  those  of  Onsager  and  Kovasnay. 

Lin  et  al^  solved  these  equations  10s  f- 

\ \ 5/ 

numerically  for  a number  of  cases.  One  \ V 3 

particular  case  is  of  special  interest  !q5  _ \ 

here,  and  it  is  given  in  Fig.  4 which  is  - \ 

copied  from  their  Fig.  6.  It  shows  both  \ 

H e I04  - \ 

the  KE  and  temperatures  spectra.  This  . \ 

example  has  the  following  values  for  the  \ 

parameters:  C4  = 0.3,  C2  1 1,  F^  - 0.01,  \ 

r 0.001,  and  F^.  0.001.  From  these  \ 

parameter  values,  it  can  be  seen  that  10  - \ 

shear  effects  in  this  example  are  less  \ 

important  than  bouyancy  effects.  The  O 10'  - \ 

figure  displays  the  -5/3  slope  in  the  KE  \ 

spectrum  at  both  high  and  low  ends,  but  10°  - \ 

a -2.  75  in  the  midrange  where  bouyancy  ' \ 

dissipation  effects  become  dominant.  Iq_i  _ \ \ 

This  brings  out  two  important  points.  - \ 

First,  it  clearly  illustrates  the  anoma- 

lous  -5/3  range  effect.  In  that  range, 

the  spectral  transfer  term,  X 5 ' ^ ' 2 


Figure  4.  Minus  5/3  Spectrum  With 
"Jog”  Due  to  Turbulence  Production 
(after  Lin  et  al,  1969a) 


is  essentially  constant  in  the  numerical  io«  |0'5  I0'2  I0'1  10°  10' 

simulation,  and  this  is  in  accordance  x = 

ko 

with  our  previous  discussion.  Second, 

this  figure  may  provide  an  explanation  Figure  4.  Minus  5/3  Spectrum  With 

for  some  of  the  spectral  nonlinearities  "Jog"  Due  to  Turbulence  Production 

• ,,  , , ..  , , , , . . (after  Lin  et  al,  1969a) 

seen  in  the  data.  No  doubt,  by  adjusting 

the  parameters  so  that  production  is 

relatively  dominant  over  bouyancy,  one  could  cause  the  intermediate  range  slope 
to  be  flatter  (rather  than  steeper)  than  -5/3.  Note  that  "jogs"  also  appear  in  some 
HICAT  data. 

It  should  be  emphasized  that  the  idealized  treatment  that  we  have  given  for 
the  spectra  of  K-H  turbulence  would  predict  a -5/3  slope  for  all  k between  the 
inner  and  outer  length  (assuming  a - 0,  that  is,  a perfect  balance  between  produc- 


*This  refers  now  to  velocity  fluctuations. 


tion  and  bouyancy).  * However,  in  reality,  neither  "a"  will  be  zero  nor  will  strengths 
or  production  and  bouyant  dissipation  interactions  (between  fluctuations  and  mean 
quantities)  match  perfectly.  Hence,  spectral  nonlinearities  should  hardly  come  as 
a surprise.  Perhaps  the  more  surprising  aspect  of  the  data  is  the  extent  to  which 
it  is  linear  in  the  majority  of  cases  published  In  any  case,  it  is  exactly  such 
nonlinearities  which  could  lead  the  unwary  to  an  erroneous  value  for  r. 

K.  l'l  RBI  I.KNT  DEf.-W  \\I)  (IKNKR  \TI()\  01  1 \H(.I  EDDIES 
8.1  Introduction 

Section  8 examines  two  further  aspects  of  stratospheric  turbulence  in  relation 
to  the  anomalous  -5/3  spectrum.  The  first  concerns  the  possibility  that  turbulence 
persists  for  a relatively  long  time  after  the  energy  input  has  stopped,  that  is,  after 
Ri  1.  It  can  then  be  argued  that  in  situ  measurements  of  stratospheric  turbulence 
are  much  more  likely  to  be  of  decaying  turbulence  rather  than  turbulence  with 
energy  input.  The  second  aspect  to  be  discussed  concerns  the  generation  of  large 
horizontal  eddies  by  the  mean  shear  during  the  "energy  production"  phase  of  the 
Kelvin-Helmholtz  event.  The  possibility  will  be  discussed  that  such  large  eddy- 
structures  could  play  the  role  of  an  outer  length  during  the  decay-  period.  This 
would  provide  an  additional  rationale  for  the  presence  of  -5/3  slopes  at  consider- 
ably 2arge  wavelengths. 

If  in  fact  the  -5/3  spectra  are  only  observed  well  beyond  the  time  when  there 
is  energy  input,  then  the  above  theory  (based  upon  compensation  between  production 
and  bouyant  dissipation)  would  be  replaced  by  a theory  based  on  the  idea  that  pro- 
duction and  bouyancy  in  the  layer  are  negligible. 

In  the  discussions  below,  we  shall  assume  that  radiation  effects  can  be 

4 Q 

neglected.  This  is  reasonable  since,  according  to  McClatchey  et  al,  there 
would  be  no  more  than  1°K  change  per  day  from  radiation  effects  in  the  strato- 
sphere. A second  assumption  is  that  atmospheric  billow  events  result  in  fairly 
complete  mixing  within  the  layer  so  that  during  decay,  mean  potential  temperature 
gradient  and  mean  shear  are  negligible  (Turner50).  This  assumption  is  important, 
not  only  for  present  purposes,  but  also  for  the  purpose  of  explaining  vertical 
transport  in  the  stratosphere  by  turbulence  (see  Rosenberg  and  Dew  an*).  It  would, 
therefore,  be  very  useful  to  validate  this  assumption  by  in  situ  measurements. 

-7 

*Note  that  there  is  a k dependence  at  large  values  of  k beyond  the  inner  scale. 

49.  McClatchey,  R.,  Kuhn,  W.R.,  Bailey,  P.  , and  Ellington,  R.  (1975)  The 

Natural  Stratosphere  of  1974  (TAP  Monograph  1.  Final  Report.  DOT. 

CIAP,  DOT -TS7 - 75 - 5 1 , p.  4-77, 

50.  Turner,  J.  S.  (1973)  Bouyancy  Effects  in  Fluids,  Cambridge  1'niversity  Press. 
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However,  one  can  interpret  the  radar  study  of  Browning  and  Watkins,  which  is 
mentioned  below,  as  lending  some  support  to  our  assumption. 

8.2  D*wa\ 

From  the  above  we  see  that  an  important  question  is  "How  long  does  strato- 
spheric turbulence  persist?"  In  the  initial  stages,  the  turbulence  consists  of  a 

breaking  K-II  wave  which  takes  less  than  5 min  to  turn  over  (5  min  is  about  the 

12 

period  for  a bouyancy  oscillation  cf  Woods  and  Wiley  ).  In  the  ocean,  where 
this  process  occurs  at  low  Reynolds  number,  one  would  estimate  the  turbulence 
duration  as  being  on  the  order  of  the  turnover  time  (Woods1).  Such  is  not  the 
case  in  the  atmosphere  due  to  high  Reynolds  number  there.  We  would  therefore  be 
interested  to  know  how  long  the  turbulence  persists  beyond  the  time  it  is  fed  by 
the  mean  shear  through  layer  spreading  etc.  (the  latter  being,  presumably  on  the 
order  of  5 min). 

First,  let  us  consider  experimental  observations.  Apparently,  the  only 
available  direct  measurements  of  turbulent  persistence  in  the  stratosphere  are 
those  of  the  H1CAT  program.  The  longest  observation  to  be  found  there  was 
14  min  and  14  sec.  In  contrast  to  the  stratosphere,  there  are  many  observations 
of  CAT  in  the  troposphere  which  are  of  interest  provided  we  remember  that  the 
turnover  time  there  is  about  double  that  of  the  stratosphere.  These  involve  both 
radar  and  radiosonde  measurements.  Boucher’’-  has  published  a radar  observa- 
tion of  tropospheric  CAT  which  occupies  a 500  m thick  layer  and  which  persisted 
for  1 hr.  Perhaps  the  most  useful  radar  study  in  this  connection  is  that  of 
Browning  and  Watkins5'-  who,  as  mentioned  above,  described  a typical  K-IT  billow 
sequence  as  consisting  of  20  min  of  visible  billow  structure  followed  by  a double 
layer  of  turbulence  persisting  for  2 hr.  This  observation  is  consistent  with  the 
idea  that  beyond  the  initial  20  min  the  turbulence  continues  for  about  2 hr  without 
bouyancy  damping  (mean  gradients  inside  t lie  layer  having  beer,  destroyed  b\ 
mixing).  At  the  top  and  bottom  edges  of  the  layer,  enough  entrainment  continues 
to  occur  with  sufficient  intensity  as  to  cause  (for  2 hr)  an  observable  reflection 
due  to  the  inhomogeneities  of  index  refraction  thus  generated.  If  this  interpretation 
is  correct,  one  can  expect  that  stratospheric  turbulence  occurs  for  hours  rather 
than  minutes. 


51.  Browning,  K.A.,  and  Watkins,  C.  D.  (1970)  Observations  of  clear  air  turbu- 

lence by  high  power  radar,  Nature  227 :2G0-263 . 

52.  Boucher,  R.  J.  (1973)  Mesoscale  history  of  a small  patch  of  clear  air  tur- 

bulence, J.  of  Appl,  Meteorol.  12:814-821. 

53.  Vinnichenko,  N.  K. , Pinus,  N.  Z.  , Shmeter,  S.  M. , and  Shut-,  G. \.  (1973) 

Turbulence  in  the  Free  Atmosphere,  Consultants  Bureau,  New  York. 
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The  radiosonde  observations  reported  in  Vinnichenko  et  al^  are  consistent 
with  the  idea  that  the  det  ay  of  atmospheric  K-H  turbulence  takes  hours.  They  say 
that  the  probability  that  the  turbulent  situation  remains  unchanged  at  12  to  20  km 
altitude  is  less  than  50  percent  over  the  period  1.  5 to  (>.  0 hr. 

In  the  troposphere,  persistence  can  be  due  to  factors  other  than  decay,  for 
example,  continually  increased  mean  shears  due  to  the  deformation  effects  of 
moving  fronts.  In  the  stratosphere  one  can  imagine  that  mountain  waves  could 
prolong  turbulence.  The  above  must  be  viewed  in  this  context. 

A crude  theoretical  estimate  of  decay  time  can  be  obtained  as  follows.  (Com- 
pare problem  1.2  of  Tennekes  and  I.umlev'").  l.et  us  imagine  a cubical  box  of 
lu  e L ( >ntaining  atmosphere  in  turbulent  motion.  Assuming  no  energy  source, 
the  turbulent  decays,  but,  since  I,  is  fixed,  we  can  assume  the  length  scale  is 
constant  and  equal  to  I..  We  derive  an  expression  for  the  decay  of  kinetic  energy, 
j 'y-  -,  as  follows: 


Integrating,  we  find  the  decay  time  is 


At 


31. 


) 


(77) 


where  u.  is  the  initial  turbulent  velocity  and  Uj.  is  the  final  velocity. 

In  this  crude  model  let  us  set  L 1000  m,  the  observed  I1ICAT  turbulence 
layer  thickness  mentioned  earlier.  For  u.  and  u^,  let  us  use  the  value  for  "strong 
turbulence  and  the  value  for  "weak”  turbulence,  respectively,  as  found  in 
Vinneshenko  et  al  p.  100  namely,  U;  2.  5 m/s  and  Uj-  0.  25  m/s.  In  this  case. 
At  ■ hi  " which  is  consistent  with  the  above  considerations.  Obviously,  At  depends 
greatly  on  what  we  choose  for  Uf  and  the  above  choice  is  perhaps  conservative. 

This  value  for  duration  disagrees  with  the  results  of  Badgley^  who  calculated 
decay  time  on  the  basis  of  constant  e.  In  general,  however,  e decreases  with 
time  and  this  prolongs  the  duration  of  decay.  On  the  other  hand,  the  value  cal- 
culated for  duration  depends  very  much  upon  the  final  turbulence  velocity;  thus,  a 
question  arises  concerning  what  final  velocity  should  be  considered  reasonable. 
Zimmerman  and  I.oving““  calculated  several  values  of  e from  the  I1ICAT  spectra 
which  ranged  from  2 f>2  to  24  cm^/sec^.  Converting  these  to  initial  and  final 
velocities  and  again  assuming  L 1000m,  we  find  At  - 20  min,  and  this  is  four 
times  longer  than  the  turn  over  time.  Thus,  even  by  this  conservative  estimate, 
it  is  reasonable  to  suppose  that  most  of  the  HICAT  turbulence  data  is  in  the  decay 
phase. 

54.  Badglev,  I . I.  (lriiip)  Large  scale  processes  contributing  energy  to  clear  air 
turbulence,  in  Clear  Air  Turbulence  and  Its  Detection,  edited  by  V.  H.  Pao 
and  A.  Goldhurg,  Plenum  Press,  New  York. 
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8.3  (iem-ration  of  Large  Horizontal  Eddies  l>y  Shear,  ami  Cascade 

There  is  a possibility  that,  in  the  anisotropic  shear  turbulence  which  we  are 
here  considering,  the  turbulence  cascade  starts  at  a scale  much  larger  than  the 
"outer  scale"  of  turbulent  layer  thickness.  If  such  were  the  case,  it  would  not  be 
possible  to  treat  this  situation  theoretically  with  the  help  of  present  theories  of 
turbulence.  On  the  other  hand,  if  such  an  effect  were  actually  present,  it  would 
further  explain  an  anomolous  -5/3  range  to  unexpectedly  large  scales.  In  the 
following,  we  shall  state  the  case  for  this  possibility.  We  shall  assume  that  the 
effect  would  be  seen  only  during  the  decaying  stage  of  the  turbulence.  But,  as  the 
above  discussion  indicates,  we  would  expect  that  this  stage  is  significantly  more 
probable  than  the  initial  "fed"  stage,  in  any  in  situ  observation.  This  implies 
that  the  anomalous  -5/3  slope  could  be  more  often  due  to  the  absence  of  both 
production  and  bouyant  dissipation  terms  rather  than  their  mutual  cancellation. 

In  such  cases,  the  measured  value  of  t would  become  reliable  even  in  the 
anomalous  range. 

Terms  VII  in  Eq.  (5(i)  and  VI  in  Eq.  (57)  describe  the  rotation  and  deformation 
effects  of  mean  shear  on  turbulence,  'i  he  latter  consist  of  compression  at  135° 
and  stretching  at  45°  where  0°  is  taken  a.->  downstream  and  90°  is  vertical,  and  a 
positive  value  for  shear  is  assumed.  This  effect  of  shear  is  duscussed  in  a 
number  of  places,  eq  Phillips^,  Townsend”’,  Eumley  Moffatt^,  and  I.umley 
and  Panofsky'"1.  In  the  above  discussion,  we  have  referred  to  it  as  a transfer  of 
energy  in  k-  space  due  to  mean  shear.  At  13  5*  , the  transfer  is  from  lower  to 
higher  wavenumbers  and  at  45°  it  is  from  higher  to  lower  wavenumbers.  Lumlej 
and  Panofsky"  1 showed  that  the  spherically  averaged  effect  is  a transfer  to  higher 
wavenumbers.  However,  the  geometry  ,,f  the  situation  clearly  indicates  that,  for 
horizontal  eddies,  the  transfer  is  to  larger  scale. 

As  discussed  in  Phillips'1, 5 the  generation  of  large  and  elongated  eddies  in  the 
horizontal  direction  is  not  only  theoretically  expected  but  has  been  observed 
experimentally  in  wall  turbulence.  Panofsky  and  Deland ' ’ ' ‘ have  studied  such 

55.  Phillips,  0.51.  (1909)  Shear-flow  Turbulence,  Annual  Rev.  of  Fluid  Dynamics 

Vol.  I.  edited  by  W.  H.  Sears,  and  M.  Van  Dyke,  Annual  Reviews,  Inc. , 
Palo  Alto,  California. 

56.  Townsend,  A.  A.  (1956)  The  Structure  of  Turbulent  Shear  I'lnw,  Cambridge 

University  Press. 

57.  Lumley,  J.  L.  (1964b)  Spectral  energv  budget  in  wall  turbulence,  Phys. 

Fluids  7:190-196. 

58.  Moffatt,  H.  K.  (1967)  The  interaction  of  turbulence  with  strong  wind  -.hear,  in 

Atmospheric  Turbulence  and  Radio  Wave  Propagation,  edited  by  A.  M. 
Yaglom,  and  V.  I.  Tatarsky,  Nauka,  Moscow. 

59.  Panofsky,  H.A.,  and  Deland,  R.J.  (1957)  Structure  of  turbulence  of  O'Neill, 

Nebraska  and  its  relation  to  the  structure  at  Brookhaven,  National 
Laboratory,  t pton,  Long  Island,  AFCRC-TR-58-201,  AD-146  843, 
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elongated  eddies  in  the  boundary  layer  with  co-spectral  techniques.  It  is  perhaps 
fair  to  say  that  large  persistent  horizontal  eddies  should  be  found  in  the  turbulent 
shear  layers  of  the  stratosphere  (having  been  generated  during  the  strong  inter- 
action, early  time  of  the  turbulence)  and  that  for  this  reason,  significant  amounts 
of  energy  may  be  found  at  scales  exceeding  layer  thickness.  During  the  decay 
phase,  shear  and  temperature  gradient  will  be  absent,  and  term  I of  Eqs.  (56) 
and  (57)  would  be  the  only  source  of  energy.  The  unanswered  question  is,  can 
some  anisotropic  analogue  of  "F"  and  "Ir,j.,j."  describe  the  cascade  processes  in 
such  a way  that  an  approximate  -5  ,'i  slope  of  the  one-dimensional  horizontal 
spectra  occur  out  to  scales  exceeding  layer  thickness  during  the  decay-’  In  other 
words,  could  the  large  eddies  generated  by  the  shear  at  the  start  of  the  turbulence 
act  later,  during  decay,  as  the  prime  energy  source  of  the  turbulence,  and  could 
an  anisotropic  cascade  result  which  approximates  the  inertial  cascade  ? Unfor- 
tunately, these  questions  must  remain  unanswered  at  present.  If  the  answers  were 
affirmative,  however,  it  would  probably  also  help  one  to  explain  why  the  slope  of 
the  vertical  velocity  spectra  are  flatter  at  large  scale  than  the  horizontal  velocity 
spectra  in  the  HICAT  data. 

( OM  1,1  SIONS 

The  HICAT  data  contain  spectra  which  have  an  approximately  -5  1 3 slope  out 

to  wavelengths  orders  of  magnitude  larger  than  expected  for  an  inertial  range. 

In  order  to  explain  this  we  have  employed  treatments  based  on  those  of  I.umley-  , 
39  25 

Gisina  , and  Lin  , all  of  which  lead  to  a -5/3  slope  for  K-Il  turbulence  the  outer 
length  when  Hi  is  close  to  unity.  A key  assumption  in  some  of  these  treatments 
was  "local  inertiality. " Since  current  theory  relies  upon  spherical  averages,  it 
cannot  explain  the  large  wavelengths  for  -5  3 that  are  observed.  To  bridge  the 
gap,  we  have  invoked  the  concept  of  Gifford  which  accounts  for  such  effects  in  the 
one-dimensional  spectra  on  the  basis  of  spatial  aliasing. 

In  addition,  some  remarks  were  made  regarding  the  generation,  b\  mean 
shears,  of  large  elongated  horizontal  eddies  which  might,  during  the  decay  phase, 
act  as  the  energy  source  for  the  turbulence.  The  speculation  was  advanced  that 
an  anisotropic  cascade  effect  might  help  to  explain  extremely  long  wavelength  -5  3 
regions. 

The  practical  impact  of  these  considerations  to  the  determination  < . f t frot 
measured  spectra  was  also  noted.  In  addition,  it  can  be  conclude-!  that  one  need 
not  resort  to  wave  theory  to  explain  -5  3 slopes  at  large  wavelengths. 
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Notes  Added 


in  Proof 


1 - The  Kelvin  Helmholtz  type  of  turbulence  discussed  in  this  paper  is  intrinsically 
a non-stationary,  developing  flow  with  weak  interaction  between  turbulence  and 
the  mean  shear  (as  has  already  been  pointed  out  in  the  text)  For  this  reason  it  is, 
strictly  speaking,  not  correct  to  omit  the  contribution  to  the  kinetic  energy  cascade 
made  by  the  decay  of  the  larger  eddies  As  mentioned  in  the  text,  when  Rf 
approaches  the  value  1,  there  would  be  virtually  no  net  input  of  energy  to  the 
cascade  if  large  eddy  decay  were  ignored,  and,  Townsend  has  shown  that  | J Fluid 
mech,  3,  3*  1-372(1938)]  Rj  could  not  reach  the  value  of  unity  but  would  be  con- 
strained by  the  relation 

(1-Rf)  * f/(uw  U'). 

The  energy  balance  with  decay  can  be  described  as  follows.  Let  D be  defined 
as  the  rate  of  kinetic  energy  input  to  the  inertial  cascade  by  the  decay  of  the  largest 
eddies  and  let  P be  turbulent  production  through  mean  shear,  B the  bouvant  dissi- 
pation, F = _ , i.  e.  the  conversion  factor  from  temperature  fluc  tuation 

T T' 

dissipation  to  potential  energy  dissipation,  the  temperature  fluctuation  dissi- 
pation rate,  and  ? the  kinetic  energy  dissipation  rate.  We  then  have 

B - F f and 

r>  + P = B + e . 
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The  second  of  these  equations  can  be  derived  from  Eq.  6 by  retaining  term 
number  I.  It  should  also  be  mentioned  that  keeping  this  term  would  modify  our 
interpretation  of  Fig.  3,  p.  30  in  the  sense  that  the  large  eddies  would  be  con- 
sidered to  be  slowly  emptying  reservoirs  of  energy  which  then  cascades  down  the 
scale,  i.  e.  down  the  more  narrow  parts  of  the  pipes. 

As  Kj  approaches  unity,  P approaches  -Ff^,  hence  e approaches  D.  We 
therefore  can  have  quasi-inertiality  at  large  wavelengths  when  D is  significant. 
Thus,  when  P and  B are  significant  in  size,  and  when  Rj  is  near  unity,  then  the 
condition  for  a -5/3  sloped  spectrum  is  de/dk  <<  D/k.  Eq.  61  for  the  spectrum 
should  be  viewed  in  this  light. 

Note  also  that  we  have  assumed  that  R^  ~ Rj.  Here  again  D would  have  to  be 
significant  in  value,  since  otherwise,  as  Towmsend  (op.  cit.  ) has  shown 

1 

— R.  = (1  - Rf)  (Rf>  ( ™ )2 

e 8 V 6w  / 

From  all  these  considerations  we  now  see  that  the  assumption  of  weak  shear- 
interaction  has  a double  significance:  (a)  the  decay  rate,  D,  must  be  assumed 
significantly  large  in  order  to  allow  quasi-inertiality,  and  (b)  a strong  interaction 

- 1 r:  / o 

would  give  a k spectrum  rather  than  a k spectrum  as  was  shown  by  Gisina, 
op.  cit. 

In  Eq.  71,  when  a = 0 we  would  have  e = 0 if  D were  equal  to  zero.  This  again 
emphasizes  the  need  for  the  condition  that  D =/=  0 if  the  -5/3  formalism  is  to  be 
consistent.  However,  we  have  not  taken  into  formal  consideration  the  input  of 
energy  due  to  layer  spreading.  This  could  be  done  superficially  by  incorporating 
that  input  in  the  term  D above.  On  the  other  hand  there  exists  no  theory  at  present 
that  gives  the  spectral  distribution  of  this  input  due  to  layer  spreading. 

2 - It  should  be  emphasized  that  the  more  sophisticated  model  of  Lin,  which  was 
^ described  in  the  text,  makes  possible  a number  of  nonlinear  spectral  shapes  which 

' would  make  a measurement  of  e,  by  means  of  the  spectrum,  incorrect.  On  the 

other  hand,  when  both  the  turbulent  production  and  bouyancy  dissipation  terms 
are  very  small,  then  such  measurements  of  f should  be  correspondingly  accurate. 
As  seen  from  the  arguments  given,  such  a state  of  affairs  may  well  be  the  more 
probable  one  for  stratospheric  measurements.  Nevertheless  the  most  feasible  way 
to  check  this  is  by  measurements  of  the  spectrum  at  sufficiently  high  k such  that 
the  wavelengths  are  well  within  the  true  inertial  range. 
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3 - The  text  mentions  the  need  to  find  experimental  evidence  for  adiabatic  lapse 
rates  in  the  stratosphere  following  turbulent  breakdown.  Fig.  6.  33  on  P.  R-7fi 
from  Loving  (1975),  [ The  Natural  Stratosphere  of  1974,  CLAP  Monograph  I, 
Final  Report,  DOT,  CLAP.,  DOT- TST-75-5 1 ] seems  to  provide  a published 
example  of  this. 


